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The aim of the ongoing relativistic heavy-ion collision 
experiments is to explore the possible plasma phase of 
QCD, the so called quark-gluon plasma (QGP). High en- 
ergy partons produced in initial partonic sub-processes in 
collisions between two heavy nuclei will lose their energy 
while propagating through the dense matter formed after 
such collisions, resulting in jet quenching. The amount 
of quenching depends upon the state of the fireball pro- 
duced and the resulting quenching pattern may be used 
for identifying and investigating the plasma phase 
In order to quantitatively understand medium modifica- 
tions of hard parton characteristics in the final state, the 
energy loss of partons in the QGP has to be determined. 
There are two contributions to the energy loss of a parton 
in a QGP: one is caused by elastic collisions between the 
partons and the other is caused by radiative losses. It is 
generally believed that the radiative loss due to multiple 
gluon radiation (see [1, H] for a review) dominates over 
the coUisional one in the ultra-relativistic case. However, 
it has been shown recently that for realistic values of the 
parameters relevant for heavy-ion collisions, there is a 
wide range of parton energies in which the magnitude of 
the coUisional loss is comparable to the radiative loss for 
heavy [1, d, d, 0| as well as for light [1, Q quark flavors. 

Earlier estimates of the coUisional energy loss in the 
QGP [i,[li were obtained by treating the medium in an 
average manner, i.e., microscopic fluctuations were ne- 
glected. However, the QGP, being a statistical system, 
is characterized by omnipresent stochastic fluctuations. 
These fluctuations couple to external perturbations and 
the response of the medium to these perturbations can 
be expressed through suitable correlation functions of 
the microscopically fluctuating variables. Now, it is well 
known that the motion of charged particles in such an en- 
vironment is stochastic in nature and resembles Brownian 
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motion. The Fokker-Planck equation provides a natural 
basis for a differential characterization of such a stochas- 
tic motion. In a homogeneous and isotropic plasma, the 
Fokker-Planck equation can be recast in the form 
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where / (p, t) is the phase-space distribution of the 
test particle. Di in ^ is the drag coefficient - a 
quantity closely related to energy loss per unit length 
dE/dx [i, 0, [III, [H, [11. It is to be noted that starting 
from 71 = 2, all terms in ([1]) are statistical in origin, i.e, 
they arise due to microscopic field fluctuations, to lead- 
ing order in [ij, [l^ . For the drag coefficient or energy 
loss, the fluctuating field as well as the polarization field 
contribute as will be seen later. The effect of field fluc- 
tuations on the passage of a charged particle through a 
non-relativistic classical plasma has been worked out by 
several authors [H, [H, 111, [13, [H, [H, [H in the past. 
This effect leads to an energy gain of the particles and 
is most effective in the low velocity limit. Given the fact 
that the subject of the energy loss is of topical inter- 
est, it is the principal motivation of the present article 
to quantitatively estimate the effect of microscopic elec- 
tromagnetic fluctuations on the energy loss of a heavy 
quark passing through an equilibrium, weakly-coupled 
QGP within the semiclassical approximation. 

The semiclassical approach was adopted earlier to cal- 
culate the coUisional energy loss of a heavy quark due to 
polarization effects of the medium [l^ . It is assumed 
that the energy lost by the particle per unit time is small 
compared to the energy of the particle itself so that the 
change in the velocity of the particle during the motion 
may be neglected, i.e, the particle moves in a straight line 
trajectory. The energy loss of a particle is determined by 
the work of the retarding forces acting on the particle in 
the plasma from the chromo-electric field generated by 
the particle itself while moving. So the energy loss of the 
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particle per unit time is given by 
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where the field is taken at the location of the particle. 
Within linear response theory the correlation function of 
the fluctuations of charge and current densities and the 
electromagnetic fields in the medium with space-time dis- 
persion are completely determined in terms of the dielec- 
tric tensor of the medium. In the Abelian approximation, 
the total chromo-electric field induced in the QGP can 
be related to the external current of the test charge by 
solving Maxwell's equations and the equation of continu- 
ity 
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with the color charge current j". It should be noted that 
in this approximation, the dielectric tensor is diagonal 
in the color indices [2l| . In an isotropic and homogeneous 
plasma the dielectric tensor etj can be decomposed into 
longitudinal and transverse parts, 



(w, k) = EL {iu, k) rt. + CT (w, k) V, 
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The 



where, Vlj = kikj/k^ and Vjj = dij — kikj/k"^ 
gauge invariant high-temperature expression for the di- 
electric functions which is in accordance with the Abelian 
approximation (see below) are given by (see e.g. [Ml, HI]) 
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where m|j = g^T^ (1 -I- Nf/6) is the Debye mass squared. 
Substitution of ^ together with (g]) and (|5l) in Q gives 
the polarization loss of the moving parton [9|, [lOj . 

The previous formula for the energy loss in ([2|) does not 
take into account the field fluctuation in the plasma and 
the particle recoil in collisions. To accommodate these 
effects it is necessary to replace ^ with [l^ [l^l , 
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where (■ ■ • ) denotes the statistical averaging operation. 
It is to be noted that here two averaging procedures are 
performed: I) an ensemble average w.r.t the equilibrium 
density matrix and II) a time average over random fluc- 
tuations in plasma. These two operations are commuting 
and only after both of them are performed the average 
quantity takes up a smooth value [3. In the following, 
we will explicitly denote the ensemble average by {■ ■ ■) ^ 
wherever required to avoid possible confusion. 



The electric field £t in ([6]) consists of the induced field £ 
given by and a spontaneously generated microscopic 

field the latter being a random function of position 
and time. 

The classical equation of motion of the particles in the 
electromagnetic field have the form, 
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Integrating the equation of motion ([7|) we find. 
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where J-t = £t + v x Bt and Eq is the initial parton 
energy. The magnetic field Bt consists of an induced 
and a microscopic part as in the case of £t. The mean 
change in the energy of the particle per unit time is given 
by ([6]). The stochastic time dependence as embodied in 
Q comes because of explicit fluctuations in time and 
motion of the particle from one field point to another. 
The dependence on the latter will be expanded about 
the mean rectilinear motion where the former will be left 
untouched. This is done by picking a time interval At 
sufficiently large with respect to the time scale of random 
electromagnetic fiuctuations in the plasma (T2) but small 
compared with the time during which the particle motion 
changes appreciably (ti). 



Ti > At > Ta 
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Parametrically, ri ~ [g"'nn(l/g)] and T2 ~ 

{gT) ^ [53. In the weak coupling limit, these scales are 
widely separated. Keeping only the leading order terms 
in the expansion we get (see eq. 8.31 in [l6l | and eq. 
13.1.3.4 in [13) 
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where, rb (0 = "^ot- Substituting (fTO)l in ([6] 
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Since the mean value of the fluctuating part of the field 

equals zero, (E^ — 0, l^£t {rit) , t)^ equals the chromo- 

electric field produced by the particle itself in the plasma. 
The first term in ((TT|) therefore corresponds to the usual 
polarization loss of the parton calculated in Q . Provided 
there exists a hierarchy of scales ([9]), it can be shown 
that the polarization field does not contribute to lead- 
ing order in the correlations functions appearing in the 
second and third terms in as in the case of higher 
order Fokker- Planck coefficients [H, [l^l ■ These terms 
correspond to the statistical change in the energy of the 
moving parton in the plasma due to the fluctuations of 
the chromo-electromagnetic fields as well as the velocity 
of the particle under the influence of this field. The sec- 
ond term in (fTTj) corresponds to the statistical part of the 
dynamic friction due to the space-time correlation in the 
fluctuations in the electrical field whereas the third one 
corresponds to the average change in the energy of the 
moving particle due to the correlation between the fluctu- 
ation in the velocity of the particle and the fluctuation in 
the electrical field in the plasma. The temporal averaging 
in pip , by definition, includes many random fluctuations 
over the mean motion. However, the correlation function 
of these fluctuations are suppressed beyond their char- 
acteristic time scales. This allows us to formally extend 
the upper limits in time-integrations in pip to infinity. 

Now within linear response theory, the power spectrum 
of the chromo-electromagnetic fields follows from the 
fluctuation-dissipation theorem and is completel y d eter- 
mincd by the dielectric functions of the medium 
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where, 
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with 77 = k/u). 

Using the Fourier transform of £i together with ([1 
and psp , we obtain from pip the energy loss of the par 
ton due to fluctuations as (see appendix). 
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where (^^'l^ a-nd (^^rj denote the longitudinal and trans- 
verse field fluctuations, respectively, and E ^ Eq is the 
initial energy of the parton. 

Eq. p4p can be recast as (see appendix). 
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where F{uj,k) — 87rtj^ImeL/|eLp and G{uj,k) = 
167rImeT/|eT ~ fc^/w^p and vq = v. This result is ob- 
viously gauge invariant if we use there the semiclassical, 
gauge invariant expression for the dielectric functions ([5]). 

The above expression gives the mean energy (per unit 
time) absorbed by a propagating particle from the heat 
bath. Physically, this arises from gluon absorption. 
Thermal absorption of gluons was also shown to reduce 
the radiative energy loss [2^. We arrive at a somewhat 
similar conclusion as there, albeit in a different context. 
It is to be noted here that since the spectral density 

of field fluctuations (^Ff^^j^ are positive for positive fre- 
quencies by definition, according to (fT5|) the particle en- 
ergy will grow due to interactions with the fluctuating 
flelds. 
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FIG. 1: (Color online) Various contributions to the energy 
loss of charm quark in the QGP. The collisional energy loss 
(dashed line) is taken from Ref.[^. 



The contribution from field fluctuations to the heavy 
quark energy loss is shown in Fig. [T] and Fig. [21 Our 
choice of the parameters is Nj = 2, T = 250 MeV, as = 
0.3, and for the charm and bottom quark masses we take 
1.25 and 4.2 GeV, respectively. For the upper integration 
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where g ~ T is the typical momentum of the thermal 
partons of the QGP. 

In Fig. [3] and Fig. 2] we show the relative collisional en- 
ergy loss of a charm and bottom quark where the effect 
of field fluctuations is taken into account. It is evident 
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FIG. 2: (Color online) Same as Fig.[T]but for a bottom quark. 
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FIG. 4: (Color online) Same as Fig.[3]but for a bottom quark. 
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FIG. 3; (Color online) Relative importance of the fluctuation 
loss compared to the collisional energy loss of Ref.[^. We take 
the path length to be L = 5 fm. 



that the effect of the fluctuations on the heavy quark en- 
ergy loss is significant at low momenta. For momenta 
4 — 20 GeV the fluctuation effect reduces the collisional 
loss by 17 — 39% for charm and 12 — 31% for bottom. At 
higher momenta, as it will be relevant for LHC, the rela- 
tive importance of the fluctuation gain to the collisional 
loss decreases gradually. The fluctuation gain diverges 
linearly as ?; ^ 0. This can be understood by looking at 
the leading log part, 
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where kmin = rriD is the Debye mass. The divergence 
here is kinematic as 4f In flnite for w = and therefore 



lilfl = if Ifl diverges for ^;^0. 

Let us summarize here the general assumptions made 
in this investigation: 

• We consider here only the collisional energy loss 
due to elastic scattering. 

• We consider here only heavy quarks (charm, bot- 
tom) which became of particular interest in recent 
experiments at RHIC @. 

• We use the semiclassical approximation which has 
been used to calculate the mean collisional energy 
loss 0. The semiclassical approximation has been 
shown to be equivalent to the Hard Thermal Loop 
approximation which is based on the weak coupling 
limit g — > 0. It allows a systematic calculation of 
the collisional energy loss [131 and is valid in the 
high-temperature limit of QCD. It also corresponds 
to neglecting the non-Abelian terms in the QCD 
equations of motion (see e.g. [1^). Hence we use 
systematically the Abelian approximation through- 
out this work. 

• We assume a constant momentum and tempera- 
ture independent coupling constant. Recently the 
collisional energy loss has been reconsidered using 
a model to include a running temperature depen- 
dent coupling constant. Although this leads to a 
different functional dependence on the parameters, 
the energy loss for realistic situations is somewhat 
larger but of similar size as in the case of a constant 
coupling [28l |. 

• We assume an equilibrated, isotropic and homoge- 
neous QGP (see below). 

• We assume an infinitely extended QGP. Recently it 
has been argued that finite size effects are negligible 
for the collisional energy loss [1^, . 
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Let us note here that the assumption of an equihb- 
rium condition necessary imphes isotropization in mo- 
mentum space. On the other hand, matter created in 
non-central heavy-ion colhsions is anisotropic to start 
with and the strong longitudinal expansion afterwards, 
at its own, brings an anisotropy into the system [30l |. The 
characteristic feature of such anisotropic systems is the 
presence of a Weibcl-typc instability [3l|, [32, HI] . It has 
been argued that assuming a turbulent, weakly coupled 
anisotropic QGP may provide a natural explanation for 
the observed rapid isotropization time |34|. for the small 
shear viscosity to entropy density ratio (35[, or for di- 
hadron correlation functions [s^ . 

The analysis presented above can in principle be ex- 
tended to the case of a non-equilibrium, anisotropic 
QGP, provided the power spectrum of the electromag- 
netic fluctuations are known. One possible way is to 
simulate them on the lattice [37| . Fluctuations are much 
stronger in non-equilibrium situations than in thermal 
systems [11]. The effect of field fluctuations on the pas- 
sage of a heavy quark is therefore expected to be stronger 
in an anisotropic QGP [s^. Interestingly, the energy loss 
for a heavy quark in an anisotropic QGP without taking 
field fluctuations into account is negative (corresponding 
to energy gain) at low energies similar to our case [40| . 
Furthermore, in the case of an anisotropic plasma there 
could be an "anti-Landau" damping mechanism which 
would lead to an energy loss by fluctuations [i^] . 

Recently heavy quark probes at RHIC have posed new 
challenges to the theoretical understanding of the parton 
energy loss. As shown by Wicks et al [3] recent measure- 
ment of the non-photonic single electron data cannot be 
explained by the radiative loss alone. If the coUisional 
energy loss is included, the agreement is better but not 
satisfactory. As seen it is obvious that the gain due to 
the field fluctuations will lead to reduction in coUisional 
energy loss and hence more disagreement with single elec- 
tron data. On the other hand, using the transport coef- 
ficients within pQCD energy loss calculations [4l| the el- 
liptic flow coefficient of single electrons V2 is limited only 
to 2 — 3% in semi-central Au-Au collisions, while the ex- 
perimental values [H] reach up to 10% around transverse 
electron momenta of 2 GeV/c. It is suggested recently 
in Ref. [i^ that not only the energy loss but also the 
energy gain in low momenta may be required for obtain- 
ing larger theoretical V2 values. It will be interesting to 
find out whether the inclusion of the fluctuation gain can 
shed light on V2 at low momenta. 



APPENDIX A: DERIVATION OF JTID 

Since the power spectrum is diagonal in color space 
according to (|T2)) we can pull out the color factor from 
the electric fields and perform the color sum in the sec- 
ond and third terms in (fTT|) easily. With Q^Q" = Cfols, 
where Ci? is the quadratic Casimir in fundamental rep- 



resentation, we obtain, e.g., for the second term in (jlip 
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Setting t — ti = T and provided that there exist scales ti 
and T2 discussed earlier, we can write (|A1[) as, 
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Expressing £ (fo (t) , t)'s in (|A2[) in Fourier modes, utiliz- 
ing the fact that the unperturbed orbit is a straight line 
trajectory Tq (t) = Vgi, and 
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Similarly for the third term in (|lip we find. 
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Using (|A3|) and the fact the field correlations vanish for 

uj — > ±cx), we get, 



It is to be noted that in contrast to (jA4[) here only the 
longitudinal part survives. This is due to the fact that 
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differentiation w.r.t rj in (jA6p brings down one power of 
kj which, operating on field correlators (|13p . removes the 
transverse part. Adding (|A4p and (jA7P we get p^ . 

Integrating the first term in with respect to the 
azimuthal angle yields, 



#1 = 



C 



dk k' 



+1 



P\k-vo 



(A8) 



Substituting w = kvQrj the 77 integration can be per- 
formed and we get 
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Now (^LoE'j^j can be written as ri^ {u) f (to) where ns (<^) 

is the Bose distribution and f (o-!) = /(— w). Using 
nB (w) = — [1 + hb (w)], and using u = kvo we can write 
[1) as; 
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Similarly, we obtain the transverse contribution as, 
Cpas 



#2 



dk k 
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